Functional protein dynamics in a crystal.

Proteins are molecular machines and to understand how they work, we need to understand how they move. New pump-probe time-resolved X-ray diffraction methods open up ways to initiate and observe protein motions with atomistic detail in crystals on biologically relevant timescales. However, practical limitations of these experiments demands parallel development of effective molecular dynamics approaches to accelerate progress and extract meaning. Here, we establish robust and accurate methods for simulating dynamics in protein crystals, a nontrivial process requiring careful attention to equilibration, environmental composition, and choice of force fields. With more than seven milliseconds of sampling of a single chain, we identify critical factors controlling agreement between simulation and experiments and show that simulated motions recapitulate ligand-induced conformational changes. This work enables a virtuous cycle between simulation and experiments for visualizing and understanding the basic functional motions of proteins.

Investigating direct current potentials that affect native protein conformation during trapped ion mobility spectrometry-mass spectrometry.

Trapped ion mobility spectrometry-time-of-flight mass spectrometry (TIMS-TOFMS) has emerged as a tool to study protein conformational states. In TIMS, gas-phase ions are guided across the IM stages by applying direct current (DC) potentials (D1-6), which, however, might induce changes in protein structures through collisional activation. To define conditions for native protein analysis, we evaluated the influence of these DC potentials using the metalloenzyme bovine carbonic anhydrase (BCA) as primary test compound. The variation of DC potentials did not change BCA-ion charge and heme content but affected (relative) charge-state intensities and adduct retention. Constructed extracted-ion mobilograms and corresponding collisional cross-section (CCS) profiles gave useful insights in (alterations of) protein conformational state. For BCA, the D3 and D6 potential (which are applied between the deflection transfer and funnel 1 [F1] and the accumulation exit and the start of the ramp, respectively) had most profound effects, showing multimodal CCS distributions at higher potentials indicating gradual unfolding. The other DC potentials only marginally altered the CCS profiles of BCA. To allow for more general conclusions, five additional proteins of diverse molecular weight and conformational stability were analyzed, and for the main protein charge states, CCS profiles were constructed. Principal component analysis (PCA) of the obtained data showed that D1 and D3 exhibit the highest degree of correlation with the ratio of folded and unfolded protein (F/U) as extracted from the mobilograms obtained per set D potential. The correlation of D6 with F/U and protein charge were similar, and D2, D4, and D5 showed an inverse correlation with F/U but were correlated with protein charge. Although DC boundary values for induced conformational changes appeared protein dependent, a set of DC values could be determined, which assured native analysis of most proteins.

Multilevel superposition for deciphering the conformational variability of protein ensembles.

The dynamics and variability of protein conformations are directly linked to their functions. Many comparative studies of X-ray protein structures have been conducted to elucidate the relevant conformational changes, dynamics and heterogeneity. The rapid increase in the number of experimentally determined structures has made comparison an effective tool for investigating protein structures. For example, it is now possible to compare structural ensembles formed by enzyme species, variants or the type of ligands bound to them. In this study, the author developed a multilevel model for estimating two covariance matrices that represent inter- and intra-ensemble variability in the Cartesian coordinate space. Principal component analysis using the two estimated covariance matrices identified the inter-/intra-enzyme variabilities, which seemed to be important for the enzyme functions, with the illustrative examples of cytochrome P450 family 2 enzymes and class A $\beta$-lactamases. In P450, in which each enzyme has its own active site of a distinct size, an active-site motion shared universally between the enzymes was captured as the first principal mode of the intra-enzyme covariance matrix. In this case, the method was useful for understanding the conformational variability after adjusting for the differences between enzyme sizes. The developed method is advantageous in small ensemble-size problems and hence promising for use in comparative studies on experimentally determined structures where ensemble sizes are smaller than those generated, for example, by molecular dynamics simulations.

Conformational free-energy landscapes of a Na+/Ca2+ exchanger explain its alternating-access mechanism and functional specificity.

Secondary-active transporters catalyze the movement of myriad substances across all cellular membranes, typically against opposing concentration gradients, and without consuming any ATP. To do so, these proteins employ an intriguing structural mechanism evolved to be activated only upon recognition or release of the transported species. We examine this self-regulated mechanism using a homolog of the cardiac Na+/Ca2+ exchanger as a model system. Using advanced computer simulations, we map out the complete functional cycle of this transporter, including unknown conformations that we validate against existing experimental data. Calculated free-energy landscapes reveal why this transporter functions as an antiporter rather than a symporter, why it specifically exchanges Na+ and Ca2+, and why the stoichiometry of this exchange is exactly 3:1. We also rationalize why the protein does not exchange H+ for either Ca2+ or Na+, despite being able to bind H+ and its high similarity with H+/Ca2+ exchangers. Interestingly, the nature of this transporter is not explained by its primary structural states, known as inward- and outward-open conformations; instead, the defining factor is the feasibility of conformational intermediates between those states, wherein access pathways leading to the substrate binding sites become simultaneously occluded from both sides of the membrane. This analysis offers a physically coherent, broadly transferable route to understand the emergence of function from structure among secondary-active membrane transporters.

Structural analysis of ATP bound to the F1-ATPase ß-subunit monomer by solid-state NMR- insight into the hydrolysis mechanism in F1.

ATP-hydrolysis-associated conformational change of the ß-subunit during the rotation of F1-ATPase (F1) has been discussed using cryo-electron microscopy (cryo-EM). Since it is worthwhile to further investigate the conformation of ATP at the catalytic subunit through an alternative approach, the structure of ATP bound to the F1ß-subunit monomer (ß) was analyzed by solid-state NMR. The adenosine conformation of ATP-ß was similar to that of ATP analog in F1 crystal structures. 31P chemical shift analysis showed that the Pα and Pß conformations of ATP-ß are gauche-trans and trans-trans, respectively. The triphosphate chain is more extended in ATP-ß than in ATP analog in F1 crystals. This appears to be in the state just before ATP hydrolysis. Furthermore, the ATP-ß conformation is known to be more closed than the closed form in F1 crystal structures. In view of the cryo-EM results, ATP-ß would be a model of the most closed ß-subunit with ATP ready for hydrolysis in the hydrolysis stroke of the F1 rotation.

PaCS-Toolkit: Optimized Software Utilities for Parallel Cascade Selection Molecular Dynamics (PaCS-MD) Simulations and Subsequent Analyses.

Parallel cascade selection molecular dynamics (PaCS-MD) is an enhanced conformational sampling method conducted as a "repetition of time leaps in parallel worlds", comprising cycles of multiple molecular dynamics (MD) simulations performed in parallel and selection of the initial structures of MDs for the next cycle. We developed PaCS-Toolkit, an optimized software utility enabling the use of different MD software and trajectory analysis tools to facilitate the execution of the PaCS-MD simulation and analyze the obtained trajectories, including the preparation for the subsequent construction of the Markov state model. PaCS-Toolkit is coded with Python, is compatible with various computing environments, and allows for easy customization by editing the configuration file and specifying the MD software and analysis tools to be used. We present the software design of PaCS-Toolkit and demonstrate applications of PaCS-MD variations: original targeted PaCS-MD to peptide folding; rmsdPaCS-MD to protein domain motion; and dissociation PaCS-MD to ligand dissociation from adenosine A2A receptor.

Main role of fractal-like nature of conformational space in subdiffusion in proteins.

Protein dynamics involves a myriad of mechanical movements happening at different time and space scales, which make it highly complex. One of the less understood features of protein dynamics is subdiffusivity, defined as sublinear dependence between displacement and time. Here, we use all-atoms molecular dynamics (MD) simulations to directly interrogate an already well-established theory and demonstrate that subdiffusivity arises from the fractal nature of the network of metastable conformations over which the dynamics, thought of as a diffusion process, takes place.

Intrinsically disordered regions (IDRs): A vague and confusing concept for protein function.

The term "intrinsically disordered region" (IDR) in proteins has been used in numerous publications. However, most proteins contain IDRs, the term refers to very different types of structures and functions, and many IDRs become structured upon interaction with other biomolecules. Thus, IDR is an unnecessary, vague, and ultimately confusing concept.

Terms and conditions.

Brian Plosky provides some context for a debate over the use of "intrinsically disordered" to describe regions of proteins.

Intrinsic disorder: A term to define the specific physicochemical characteristic of protein conformational heterogeneity.

In his commentary in this issue of Molecular Cell,1 Struhl reasons that the term "intrinsically disordered regions" represents a vague and confusing concept for protein function. However, the term "intrinsically disordered" highlights the important physicochemical characteristic of conformational heterogeneity. Thus, "intrinsically disordered" is the counterpart to the term "folded, " with neither term having specific functional implications.

How nuclear receptors transition between active and inactive forms: An energetic perspective.

Nuclear receptors regulate transcriptional programs in response to the binding of natural and synthetic ligands. These ligands modulate the receptor by inducing dynamic changes in the ligand binding domain that shift the C-terminal helix (H12) between active and inactive conformations. Despite decades of study, many questions persist regarding the nature of the inactive state and how ligands shift receptors between different states. Here, we use molecular dynamics (MD) simulations to investigate the timescale and energetic landscape of the conformational transition between inactive and active forms of progesterone receptor (PR) bound to a partial agonist. We observe that the microsecond timescale is insufficient to observe any transitions; only at millisecond timescales achieved via accelerated MD simulations do we find the inactive PR switches to the active state. Energetic analysis reveals that both active and inactive PR states represent energy minima separated by a barrier that can be traversed. In contrast, little or no transition is observed between active and inactive states when an agonist or antagonist is bound, confirming that ligand identity plays a key role in defining the energy landscape of nuclear receptor conformations.

Role of Protonation States in the Stability of Molecular Dynamics Simulations of High-Resolution Membrane Protein Structures.

Classical molecular dynamics (MD) simulations provide unmatched spatial and time resolution of protein structure and function. However, the accuracy of MD simulations often depends on the quality of force field parameters and the time scale of sampling. Another limitation of conventional MD simulations is that the protonation states of titratable amino acid residues remain fixed during simulations, even though protonation state changes coupled to conformational dynamics are central to protein function. Due to the uncertainty in selecting protonation states, classical MD simulations are sometimes performed with all amino acids modeled in their standard charged states at pH 7. Here, we performed and analyzed classical MD simulations on high-resolution cryo-EM structures of two large membrane proteins that transfer protons by catalyzing protonation/deprotonation reactions. In simulations performed with titratable amino acids modeled in their standard protonation (charged) states, the structure diverges far from its starting conformation. In comparison, MD simulations performed with predetermined protonation states of amino acid residues reproduce the structural conformation, protein hydration, and protein-water and protein-protein interactions of the structure much better. The results support the notion that it is crucial to perform basic protonation state calculations, especially on structures where protonation changes play an important functional role, prior to the launch of any conventional MD simulations. Furthermore, the combined approach of fast protonation state prediction and MD simulations can provide valuable information about the charge states of amino acids in the cryo-EM sample. Even though accurate prediction of protonation states in proteinaceous environments currently remains a challenge, we introduce an approach of combining pKa prediction with cryo-EM density map analysis that helps in improving not only the protonation state predictions but also the atomic modeling of density data.

Small angle X-ray scattering analysis of thermophilic cytochrome P450 CYP119 and the effects of the N-terminal histidine tag.

Combining size exclusion chromatography-small angle X-ray scattering (SEC-SAXS) and molecular dynamics (MD) analysis is a promising approach to investigate protein behavior in solution, particularly for understanding conformational changes due to substrate binding in cytochrome P450s (CYPs). This study investigates conformational changes in CYP119, a thermophilic CYP from Sulfolobus acidocaldarius that exhibits structural flexibility similar to mammalian CYPs. Although the crystal structure of ligand-free (open state) and ligand-bound (closed state) forms of CYP119 is known, the overall structure of the enzyme in solution has not been explored until now. It was found that theoretical scattering profiles from the crystal structures of CYP119 did not align with the SAXS data, but conformers from MD simulations, particularly starting from the open state (46 % of all frames), agreed well. Interestingly, a small percentage of closed-state conformers also fit the data (9 %), suggesting ligand-free CYP119 samples ligand-bound conformations. Ab initio SAXS models for N-His tagged CYP119 revealed a tail-like unfolded structure impacting protein flexibility, which was confirmed by in silico modeling. SEC-SAXS analysis of N-His CYP119 indicated pentameric structures in addition to monomers in solution, affecting the stability and activity of the enzyme. This study adds insights into the conformational dynamics of CYP119 in solution.

In silico design of misfolding resistant proteins: the role of structural similarity of a competing conformational ensemble in the optimization of frustration.

Most state-of-the-art in silico design methods fail due to misfolding of designed sequences to a conformation other than the target. Thus, a method to design misfolding resistant proteins will provide a better understanding of the misfolding phenomenon and will also increase the success rate of in silico design methods. In this work, we optimize the conformational ensemble to be selected for negative design purposes based on the similarity of the conformational ensemble to the target. Five ensembles with different degrees of similarity to the target are created and destabilized and the target is stabilized while designing sequences using mean field theory and Monte Carlo simulation methods. The results suggest that the degree of similarity of the non-native conformations to the target plays a prominent role in designing misfolding resistant protein sequences. The design procedures that destabilize the conformational ensemble with moderate similarity to the target have proven to be more promising. Incorporation of either highly similar or highly dissimilar conformations to the target conformation into the non-native ensemble to be destabilized may lead to sequences with a higher misfolding propensity. This will significantly reduce the conformational space to be considered in any protein design procedure. Interestingly, the results suggest that a sequence with higher frustration in the target structure does not necessarily lead to a misfold prone sequence. A successful design method may purposefully choose a frustrated sequence in the target conformation if that sequence is even more frustrated in the competing non-native conformations.

Switching GoÌ -Martini for Investigating Protein Conformational Transitions and Associated Protein-Lipid Interactions.

Proteins are dynamic biomolecules that can transform between different conformational states when exerting physiological functions, which is difficult to simulate using all-atom methods. Coarse-grained (CG) Go̅-like models are widely used to investigate large-scale conformational transitions, which usually adopt implicit solvent models and therefore cannot explicitly capture the interaction between proteins and surrounding molecules, such as water and lipid molecules. Here, we present a new method, named Switching Go̅-Martini, to simulate large-scale protein conformational transitions between different states, based on the switching Go̅ method and the CG Martini 3 force field. The method is straightforward and efficient, as demonstrated by the benchmarking applications for multiple protein systems, including glutamine binding protein (GlnBP), adenylate kinase (AdK), and ß2-adrenergic receptor (ß2AR). Moreover, by employing the Switching Go̅-Martini method, we can not only unveil the conformational transition from the E2Pi-PL state to E1 state of the type 4 P-type ATPase (P4-ATPase) flippase ATP8A1-CDC50 but also provide insights into the intricate details of lipid transport.

Comprehensive Analysis of Coupled Proline Cis-Trans States in Bradykinin Using ωBP-REMD Simulations.

It is well-known that proline (Pro) cis-trans isomerization plays a decisive role in the folding and stabilization of proteins. The conformational coupling between isomerization states of different Pro residues in proteins during conformational adaptation processes is not well understood. In the present work, we investigate the coupled cis-trans isomerization of three Pro residues using bradykinin (BK), a partially unstructured nonapeptide hormone, as a model system. We use a recently developed enhanced-sampling molecular dynamics method (ω-bias potential replica exchange molecular dynamics; ωBP-REMD) that allows us to exhaustively sample all combinations of Pro isomer states and obtain converged probability densities of all eight state combinations within 885 ns ωBP-REMD simulations. In agreement with experiment, the all-trans state is seen to be the preferred isomer of zwitterionic aqueous BK. In about a third of its structures, this state presents the characteristic C-terminal ß-turn conformation; however, other isomer combinations also contribute significantly to the structural ensemble. Unbiased probabilities can be projected onto the peptide bond dihedral angles of the three Pro residues. This unveils the interdependence of the individual Pro isomerization states, i.e., a possible coupling of the different Pro isomers. The cis/trans equilibrium of a Pro residue can change by up to 2.5 kcal·mol-1, depending on the isomerization state of other Pro residues. For example, for Pro7, the simulations indicate that its cis state becomes favored compared to its trans state when Pro2 is switched from the trans state to the cis state. Our findings demonstrate the efficiency of the ωBP-REMD methodology and suggest that the coupling of Pro isomerization states may play an even more decisive role in larger folded proteins subject to more conformational restraints.

Hairpin trimer transition state of amyloid fibril.

Protein fibril self-assembly is a universal transition implicated in neurodegenerative diseases. Although fibril structure/growth are well characterized, fibril nucleation is poorly understood. Here, we use a computational-experimental approach to resolve fibril nucleation. We show that monomer hairpin content quantified from molecular dynamics simulations is predictive of experimental fibril formation kinetics across a tau motif mutant library. Hairpin trimers are predicted to be fibril transition states; one hairpin spontaneously converts into the cross-beta conformation, templating subsequent fibril growth. We designed a disulfide-linked dimer mimicking the transition state that catalyzes fibril formation, measured by ThT fluorescence and TEM, of wild-type motif - which does not normally fibrillize. A dimer compatible with extended conformations but not the transition-state fails to nucleate fibril at any concentration. Tau repeat domain simulations show how long-range interactions sequester this motif in a mutation-dependent manner. This work implies that different fibril morphologies could arise from disease-dependent hairpin seeding from different loci.

In Vitro and In Silico Explorations of the Protein Conformational Changes of Corynebacterium glutamicum MshA, a Model Retaining GT-B Glycosyltransferase.

MshA is a GT-B glycosyltransferase catalyzing the first step in the biosynthesis of mycothiol. While many GT-B enzymes undergo an open-to-closed transition, MshA is unique because its 97° rotation is beyond the usual range of 10-25°. Molecular dynamics (MD) simulations were carried out for MshA in both ligand bound and unbound states to investigate the effect of ligand binding on localized protein dynamics and its conformational free energy landscape. Simulations showed that both the unliganded "opened" and liganded "closed" forms of the enzyme sample a wide degree of dihedral angles and interdomain distances with relatively low overlapping populations. Calculation of the free energy surface using replica exchange MD for the apo "opened" and an artificial generated apo "closed" structure revealed overlaps in the geometries sampled, allowing calculation of a barrier of 2 kcal/mol for the open-to-closed transition in the absence of ligands. MD simulations of fully liganded MshA revealed a smaller sampling of the dihedral angles. The localized protein fluctuation changes suggest that UDP-GlcNAc binding activates the motions of loops in the 1-l-myo-inositol-1-phosphate (I1P)-binding site despite little change in the interactions with UDP-GlcNAc. Circular dichroism, intrinsic fluorescence spectroscopy, and mutagenesis studies were used to confirm the ligand-induced structural changes in MshA. The results support a proposed mechanism where UDP-GlcNAc binds with rigid interactions to the C-terminal domain of MshA and activates flexible loops in the N-terminal domain for binding and positioning of I1P. This model can be used for future structure-based drug development of inhibitors of the mycothiol biosynthetic pathway.

N-Glycosylation as a Modulator of Protein Conformation and Assembly in Disease.

Glycosylation, a prevalent post-translational modification, plays a pivotal role in regulating intricate cellular processes by covalently attaching glycans to macromolecules. Dysregulated glycosylation is linked to a spectrum of diseases, encompassing cancer, neurodegenerative disorders, congenital disorders, infections, and inflammation. This review delves into the intricate interplay between glycosylation and protein conformation, with a specific focus on the profound impact of N-glycans on the selection of distinct protein conformations characterized by distinct interactomes-namely, protein assemblies-under normal and pathological conditions across various diseases. We begin by examining the spike protein of the SARS virus, illustrating how N-glycans regulate the infectivity of pathogenic agents. Subsequently, we utilize the prion protein and the chaperone glucose-regulated protein 94 as examples, exploring instances where N-glycosylation transforms physiological protein structures into disease-associated forms. Unraveling these connections provides valuable insights into potential therapeutic avenues and a deeper comprehension of the molecular intricacies that underlie disease conditions. This exploration of glycosylation's influence on protein conformation effectively bridges the gap between the glycome and disease, offering a comprehensive perspective on the therapeutic implications of targeting conformational mutants and their pathologic assemblies in various diseases. The goal is to unravel the nuances of these post-translational modifications, shedding light on how they contribute to the intricate interplay between protein conformation, assembly, and disease.

Mechanistic Insights Revealed by YbtPQ in the Occluded State.

Recently, several ATP-binding cassette (ABC) importers have been found to adopt the typical fold of type IV ABC exporters. Presumably, these importers would function under the transport scheme of "alternating access" like those exporters, cycling through inward-open, occluded, and outward-open conformations. Understanding how the exporter-like importers move substrates in the opposite direction requires structural studies on all the major conformations. To shed light on this, here we report the structure of yersiniabactin importer YbtPQ from uropathogenic Escherichia coli in the occluded conformation trapped by ADP-vanadate (ADP-Vi) at a 3.1 Å resolution determined by cryo-electron microscopy. The structure shows unusual local rearrangements in multiple helices and loops in its transmembrane domains (TMDs). In addition, the dimerization of the nucleotide-binding domains (NBDs) promoted by the vanadate trapping is highlighted by the "screwdriver" action at one of the two hinge points. These structural observations are rare and thus provide valuable information to understand the structural plasticity of the exporter-like ABC importers.